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We report on first principles calculations of superconductivity in a single layer of lead on a silicon 
substrate including a full treatment of phononic and RPA screened coulomb interactions within the 
parameter free framework of Density Functional Theory for superconductors. A thorough inves¬ 
tigation shows that several approximations that are commonly valid in bulk systems fail in this 
constrained 2D geometry. The calculated critical temperature turns out to be much higher than the 
experimental value of 1.86K. We argue that the only plausible explanation for the experimental T c 
suppression is the onset of fluctuations of the superconducting order parameter. 


Nature shows a clear correlation between superconduc¬ 
tivity and dimensionality as all superconductors with a 
high critical temperature (T c ), cuprates, pnictides and 
MgB 2 , have sharp two-dimensional properties. Under¬ 
standing this connection is among the most important 
targets in contemporary solid state research. It is likely 
that relevant physical mechanisms work differently in re¬ 
duced dimensionality and that approximations and theo¬ 
retical methods developed through the experience accu¬ 
mulated on three-dimensional systems have to be modi¬ 
fied for constrained geometries. Moreover fluctuation in¬ 
stabilities of the order parameter may play an important 
rolep]. 

The conclusive test to check the theoretical under¬ 
standing is to perform ab initio calculations and compare 
directly with experiments. Currently, such a test cannot 
be done for pnictides and cuprates as the pairing mech¬ 
anism is still under debate. However, it can be done for 
phononic superconductors. 

In this work we present the results of this type of anal¬ 
ysis for lead, as this phononic superconductor is correctly 
described in bulk by ab-initio methods 013] and is exper¬ 
imentally realized in the two-dimensional limit by depo¬ 
sition on a silicon 111 substrate[4}£7j. 

This Si-Pb system is constructed as shown in Fig. |TJ 
We model the Si substrate by a 111 oriented slab, which 
is passivated on the opposite side of the lead surface using 
hydrogen [8] |9 . A relatively large width of five Si-bilayers 
is chosen in order to reduce spurious size effects of the 
substrate on the Pb layer. For the same reason we con¬ 
strain the hexagonal (xy) Si unit cell to its bulk size. 
Lead is placed in the so-called striped incommensurate 
(SIC) configuration. Since we work with periodic bound¬ 
ary conditions, a vacuum of ^8 A separates the periodic 
replica of the system. Within these constraints a full re¬ 
laxation is performed. Relaxations, electronic structure, 
phonons and electron-phonon interactions have been cal¬ 
culated within Kohn-Sham (KS) density functional the¬ 
ory (DFT) [10] . The calculated electron phonon coupling 
strength results in A = 0.78. If we use the McMil¬ 
lan formula|T4] with a standard value for the parameter 
H* = 0.10 we obtain an estimation for the critical temper¬ 


ature of 1.98K. This is in very good agreement [15] with 
the experimental T c of 1.86K[4]. With this result one has 
to conclude that superconductivity in this 2D limit can 
be understood from the electronic coupling alone and no 
fluctuations are necessary to explain the physics of this 
system. 



Figure 1. (color online) SIC configuration of Pb on the Si 
(111) substrate. On the left we present the top view, and on 
the right the side view. Black lines mark the simulated unit 
cell. 


Is this really the end of the story? Is it correct to as¬ 
sume the validity approximations known to work well in 
bulk superconductors also for this low dimensional sys¬ 
tem? The answer is no, and to show this we proceed to 
deeper investigation. 

In order to avoid any adjustable parameter (as the 
above named /i*) we use density functional theory for 
superconductors [2] H6U25] (SCDFT), where electronic 
and phononic couplings are included on the same footing. 

Electronic and phononic properties - A very relevant 
property in the electronic structure (Fig. [2]) is the pres¬ 
ence of both Pb and substrate metallic bands. This 
means that Pb deposition acts as a dopand to the Si 
substrate which develops a surface metallic region. This 
metallic region fades away within a few layers. The 
presence of this additional metallic band is relevant for 
two reasons. First it may provide a contribution to 
the electron-phonon coupling and, second, it may sta¬ 
bilize fluctuations of the order parameter of the su¬ 
perconducting phase by effectively enhancing the three- 
dimensionality of the condensate. These Si metallic 
bands can be removed by using an n-doped substrate. 
We explicitly consider this case by substituting one Si 
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x pb ’ pb 

^Pb,Si 

^ Si,Si 

Xav 

max[Ai] Np b (0) Nsi(0) 


T c A p 6 (0) A Si (0) 

7c 

undoped 

0.95 

0.13 

0.06 

0.78 

0.98 0.97 0.60 

undoped 

3.42 0.71 0.32 

2.01 

doped 

1.03 

0.00 

0.00 

1.03 

1.03 1.07 0.00 

doped 

3.54 0.74 

2.74 


Table I. Electron phonon coupling coefficients. X l,j is Table II. Calculated critical temperatures, T c (in K), within 

the Fermi Surface sheet resolved coupling matrix. X av = SCDFT and superconducting gap, A (in meV), on the Lead 

iv^o) Ni(0) is the average electron phonon coupling and Si Fermi surfaces. T* is the critical temperature esti- 

where W(0) are the Fermi surface resolved DOS and N( 0) is mated using an average coupling on the Fermi surfaces, ig- 

the total DOS. max[A*] is the maximum eigenvalue of A that noring the energy dependence of dos and screened coulomb 

in BCS acts as the effective pairing to determine the critical interactions (corresponding to a g* like approximation), 

temperature [26] . 


atom (in the deep bulk) with a virtual mixture of P and 
Si, corresponding to a doping of 1 part per 240 Si. Dop¬ 
ing has a small effect on the filling level of the Pb bands, 
but completely saturates the Si- hole pockets (see Fig.j2|. 
This doped system is experimentally realized [4] and al¬ 
lows for direct comparison with results obtained in this 
study. 
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Figure 2. (color online) Band Structure (left) and density of 
states (right) near the Fermi Energy. Thick lines correspond 
to a phosphorus doped system (1 part per 224), and the color 
scale corresponds to the projection on the Pb states. Dashed 
lines are the bands in the undoped system. The P doping has 
a negligible effect on Pb bands while it completely fills the 
silicon hole pockets (indicated by green arrows). 

The degree of two-dimensionality of the SC phase, i.e. 
how the condensate extends into the substrate, is deter¬ 
mined by the lead-substrate interaction. We can distin¬ 
guish three main effects that describe how the Pb surface 
and the substrate are coupled: chemical hybridization, 
electron phonon coupling and Coulomb interaction. 

The chemical hybridization between surface and sub¬ 
strate states can be made visible by projecting the KS 
states on the Pb atomic orbitals. This analysis shows 
that the KS states near the Fermi energy are either lo¬ 
cated in the lead surface or inside the silicon bulk, with 
no overlap (see Fig. |2j. 

The electron-phonon coupling is computed for the KS 
system via linear response [10 j. Phonons may generate 


pairing between bulk and surface states. In Tab. [T] we re¬ 
port the FS-resolved el-ph coupling[27j. By considering 
the average coupling and ignoring the energy dependence 
of density of states and screened coulomb interactions (by 
approximating them with the value at the Fermi energy) 
we have a formal equivalence of SCDFT with the McMil¬ 
lan method. The resulting critical temperature of lead 
on the undoped substrate is T c = 2.OIK and T c rises to 
2.74K for the doped Si substrate. The difference in crit¬ 
ical temperatures between the doped and undoped sys¬ 
tem is caused by the fact that the undoped material has 
an mean coupling which is much weaker than the lead- 
lead intra-surface coupling alone. This implies that the 
isotropic approximation is unjustified and leads to an un¬ 
derestimation of T c . Multiband-superconductivity must 
be explicitly accounted for as in the well known case of 
MgB 2 . 

Moreover the electrons are subject to a screened 
Coulomb scattering which we treat within the RPA[28| 
[29]. This kind of interaction in bulk materials is often 
overlooked, since, acting both as a repulsive (directly) 
and attractive interaction (via Coulomb renormalization 
mechanism [30H33] ) it appears very often to be largely 
material independent. This shows up in Eliashberg based 
methods [HI j34] in the well-known rule of thumb to take 
pb* ~ 0.1. A crucial advantage of SCDFT is that via 
the matrix elements of the RPA-screened Coulomb inter¬ 
action the Coulomb renormalization effect is explicitly 
calculated, making the use of empirical parameters like 
/i* obsolete. A metallic layer on a semiconducting sub¬ 
strate is conceptually different from a bulk in that, due to 
the lower dimensionality, there is a reduced phase space 
for low energy Coulomb scattering, that is repulsive for 
Cooper pairing (in s-wave), while the space for high en¬ 
ergy scattering is not restricted, owing to the presence of 
the substrate. Therefore the Coulomb renormalization is 
unusually large in this type of system. 

Discussion - The computed critical temperature for 
the undoped (doped) system as given in Tab. [IT] is 
3.42(3.54)K. 

We have then shown that, releasing several unjusti¬ 
fied approximations, the estimated critical temperature 
of 3.54K (doped system) is far too high as compared to 
the values experimentally observed 1.86K[4], 1.5K|6] and 
1.1K0. 
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What is the source of this mismatch? To answer this 
question we have to carefully investigate the effects not 
considered in the above analysis and their possible influ¬ 
ence on superconductivity. I) We have assumed the RPA 
represent the screened Coulomb interaction. This is re¬ 
liable in the high-density limit when screening is good. 
Therefore the Pb layer is expected to be well described. 
The approximation may be less accurate for the silicon 
hole band, since these states have a low density and, thus, 
will be poorly screened. However, the strong Coulomb 
repulsion will prevent a significant contribution to su¬ 
perconductivity, therefore this inaccuracy cannot affect 
the estimated T c significantly. Surely not for the doped 
system where these bands do not even cross the Fermi 
level. II) In general, when computing the electron phonon 
pairing, vertex corrections can be safely dropped, due 
to Migdal’s theorem 34'. [35]. The shape of the Si hole 
pocket band might invalidate this conclusion. However, 
this cannot have a significant influence on the calculation 
of superconductivity in this system since, as discussed 
above, this band effectively does not take part in the 
condensation. Migdal’s theorem is also not applicable in 
the small q limit. This does not affect the estimation 
of the phononic pairing, due to the small fraction of the 
Brillouin zone in which the problem occurs. Neverthe¬ 
less we have to keep in mind that the low q physics are 
not correctly described under this assumption. Ill) In 
our calculations we do not include spin-orbit coupling ef¬ 
fects. These have been shown to be relevant both for bulk 
lead [36] [37] and lead multilayers [SBj. However the effect 
systematically increases the coupling strenght, therefore 
it can not explain our overestimation of T c . And actually 
its inclusion would lead to an even higher critical tem¬ 
perature. IV) In our work we consider only a statically 
screened Coulomb interaction. The result of dynamic 
(plasmonic) effects could lead to important modifications 
of the dielectric screening in the case of low energy sur¬ 
face plasmons. However, as first pointed out by Takada, 
this effect is known to give a positive contribution to su¬ 
perconductivity (enhancement of coulomb renormaliza¬ 
tion by the plasmonic peak[23] [39] 00] ). Therefore, if 
relevant, it would lead to a higher estimate of T c . V) An¬ 
other questionable approximation is the use of the LDA 
in the low dimensional limit. This issue has been inves¬ 
tigated in detail by Pollack and Perdew[4T| showing that 
LDA performs well as soon as the ratio between the layer 
thickness and the r s coefficient of the gas is ~ 2. In our 
case this ratio can be estimated to be of the order of 5 
and we expect the LDA to perform as reliably as usual. 
VI) Due to the poor metal-substrate coupling, the calcu¬ 
lated single particle excitation spectrum of Si presents a 
fundamental gap that is about one half of the observed 
gap in bulk silicon. This may lead to an overestimation 
in the Coulomb renormalization, and then in an overesti¬ 
mation of T c . We have therefore accounted for this effect 
in our calculations by including a scissor correction on 
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Figure 3. (color online) Real-space structure of the SC or¬ 
der parameter y(R, 6) normalized to its maximal value of 
0.0002765. As it can be interpreted as the wave function of 
condensed pairs the confinement of the SC phase to the Pb 
layer is clearly visible. The dark blue in the substrate indi¬ 
cates that it takes part in the Coulomb renormalization and 
thus reflects a proximity effect that extends throughout the 
whole substance. 


the Si bands and the resulting effect on T c correction is 
< 0.1K. 

We believe that we have considered all relevant elec¬ 
tronic pairing effects. In the bulk limit the critical tem¬ 
perature in SCDFT, using the same approximations as 
for the slab is 6.3 K that compares well with the experi¬ 
mental value of 7.2 K. 

The only mechanism that is not included in our sim¬ 
ulations and that, according to model calculations, is 
strongly suspected to suppress superconductivity, is the 
onset of fluctuations in the order parameter. While this 
could be in principle captured in SCDFT, the presence 
of infrared collective excitations of the order parameter 
is not accounted for in the present functionals. Owing to 
Mermin-Wagner’s theorem|l, 42] these fluctuations com¬ 
pletely forbid superconductivity in a strictly two dimen¬ 
sional system. In 3D systems of constrained geometry 
(such as surfaces) model calculations show that these 
fluctuations may still be relevant in the limit in which the 
thickness is of the atomic scale and the in-plane dimen¬ 
sion of the system is macroscopic03]. Due to the strong 
confinement of the SC phase to the lead layer, as is clearly 
seen in the real space structure of the order parameter of 
a one would expect to be in a regime where these 
fluctuation effects of the superconducting order parame¬ 
ter are relevant. While, as mentioned, neglecting vertex 
corrections has probably little effect on the phononic el- 
el coupling, effective interactions in the superconducting 
Nambu channel in the sense of the fluctuation propaga¬ 
tor 04] can be very important. The disagreement be¬ 
tween the calculated and experimental critical temper¬ 
ature then strongly suggests that T c is experimentally 
limited by the fluctuation regime. The superconducting 
phase rapidly stabilizes with an increasing number of Pb 
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layers;IS 07] strengthening this conclusion. 

To summarize, we report a first-principles calculation 
of the superconducting ground state of a single lead 
layer deposited on a Si (111) substrate. We account for 
phonon mediated (via linear response DFT) and screened 
coulomb pairing (RPA) within the parameter free frame¬ 
work of Superconducting Density Functional Theory. We 
have shown that the isotropic approximation is not valid 
in this surface configuration, and in particular that the 
isotropic /a* approximation used for bulk superconductiv¬ 
ity leads to a large underestimation of the critical tem¬ 
perature. Our calculations predict a critical temperature 
about 80% larger than observed in experiment. Our anal¬ 
ysis strongly suggests that this mismatch is attributed to 
the onset of long wavelength phase fluctuations of the 
superconducting order parameter. 
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